Valence-band effective-potential evolution for coupled holes 
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Abstract 

We present the metamorphosis in the effective-potential profile of layered heterostructures, for 
several III-V semiconductor binary compounds, when the band mixing of light and heavy holes 
increases. A root-locus-like procedure, is directly applied to an eigenvalue quadratic problem 
obtained from a multichannel system of coupled modes, in the context of multiband effective 
mass approximation. By letting grow valence-band mixing, it is shown the standard fixed-height 
rectangular potential-energy for the scatterer distribution, to be a reliable test-run input for heavy 
holes. On the contrary, this scheme is no longer valid for light holes and a mutable effective band 
offset profile has to be considered instead, whenever the in-plane kinetic energy changes. 
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I. INTRODUCTORY OUTLINES 



For holes, the single-band effective mass approximation becomes inappropriate for de- 
scribing the quantum properties of band-mixed states. These quasi-particles are essentially 
mixed even far away from the scattering potential we are interested in. Many theoretical 
models have been developed to face this problem in multi-layered systems. In this work, 
we deal with the k • p two-band effective-mass Kohn-Liittinger (KL) model PQ, due to its 
widely accepted accuracy for describing dynamics of elementary excitations and electronic 
properties in the valence band [21 [3] . 

We have focused the multiband-hole band mixing-phenomenon itself due to the strong 
dependency of hole quantum transport physics upon the wave vector transversal to the main 
direction of transmission (k t )- This phenomenology, early quoted by Wessel and Altarelli 
in resonant tunneling [I], has been lately stressed for real-life technological devices [SJ. 
Fundamental condensed- matter studies [SI CTJ , had propelled us into the present modelling, 
since they have predicted the modification of the effective potential in the electronic case. 
Here the key lies in the fact, that the potential-energy profile distribution in either of the 
binary-alloy slab might evolve, depending on the value of the transversal component of 
the wave vector [SI [7J. With respect to the electrons quantum transport through layered 
heterostructures, the situation become more cumbersome, whenever the holes are involved. 
These last charge carriers — been the heaviest ones — , notably dictate the threshold response 
of technological devices. 

We recall earlier contributions to the hole subband structure in quantum wells, that has 
been previously reported within the envelope-function approximation for KL Hamiltonians 
[21 [3]. The hole-state quantization calculated in truncated parabolic confining potential 
[2], as well as for different grown directions [3], were undoubtedly striking results, whose 
validity was amply discussed in the literature. Having focused mainly the hole spectrum in 
quantum systems, they do not address, for example, the problem of the effective potential 
metamorphosis that undergo the mixing of holes. While these studies [21 E] have added 
substantial contributions to the elucidation of the valence-subband structure in quantum 
systems, there remain some aspects which do not appear to have received yet sufficient 
attention and/or because of their interest deserve further clarification. This is essentially 
the case of the carriers' transverse motion influence on the effective scattering potential 
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they interact with. This is crucial for quantum transport calculation, a question soon to be 
considered partially in this paper. 

On general grounds, for k t the band offset of the effective potential V e s is given by 
the difference for 3D band edge levels [3 IB] • By letting k t grow, the band mixing effects arise 
and the effective band offset changes. A comprehensive analysis of this subject, describes 
largest reduction for the piecewise constant V e s height as a function of k t , for light holes 
(Ih) respect to that for heavy holes (hh) j8]. Several new and peculiar features for each hole 
flavor have been found, and they enhance the novelty of the present study concerning the 
results reported in the past [BHE]. We will refer to root-locus-like terminology from now on 
throughout the paper, whenever we proceed to graph on the complex plane, the eigenvalues 
evolution for the non-linear problem as the band mixing parameter changes. This complex- 
plane plot procedure, when a system's parameter varied, is amply known in the literature 
as the root locus [10] . Taking advantage of the root-locus-like technique simulations |9J, we 
were able to foretell features of the particle-scatterer interaction. 

It is worthwhile to recall, there had also been several theoretical attempts to spread 
light to the effective potential dependency as k t increases [TT1 IT2] . However, we do believe 
that "monitoring" the quadratic eigenvalue problem via the root-locus-like procedure [9], 
and plotting the metamorphosis of V e g, seem a promising way to get a better insight into 
a complicated phenomenon referred to as valence-band mixing. The last, is the target of 
the present theoretical study. In the next section, we present briefly the basic theoretical 
workbench to study V e Q in the valence band. In the following section, for highly specialized 
III-V binary-compound semiconductor, the model is numerically tested for the evolution 
of V e g . We devote the last section, to draft some conclusions. 

II. BASIC THEORETICAL FORMULATION 

An effective potential, is found useful to describe valence-band mixing in the EFA, [P2] . 
In the case envisioned here, to determine the operator W c g for the effective band offset 
potential, suffices to use the Kohn-Luttinger (KL) model Hamiltonian, [TJ, considering the 
transverse quasi-momentum k t = k x e x + k y e y , because this is the direction of the Brillouin 
Zone where is described the present KL Hamiltonian. We assume understood any modi- 
fication of the selected Brillouin Zone direction, as a change in the Hamiltonian model to 
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FIG. 1. Panel (a)/(b) displays the metamorphosis of the effective potential profile V e Q for Ih/hh 
(red/blue lines), as a function of k t and layer dimension for a GaAs -AlAs-GaAs heterostructure. 

use. The system's quantal heterogeneity is considered along the z axis, taken perpendicular 
to the heterostructure interfaces. The operator VJ^ fj, is nothing but somewhat arbitrary 
convention, valid as long as one get holds of all potential-like energy terms from the original 
Hamiltonian operator, which are z-component momentum free [61, [8] - Then 

W u W 12 " 
w? 2 W 22 

^ 22 W 12 

W£j W n _ 
is suitable for going through a standard calculation 



W eS = 



w cS - v cS u 



*(z) 



o 4 



(2) 



leading us to the effective potential band offset V e s, "felt" in some sense, by holes during 
their passage trough the heterostructure, as k t changes. 



In the theoretical formulation, we use the periodic part of the Bloch functions whose 
components u (r) = (ui, U2, 1*3, U4) 



M i = 114) = ^\( x + i y) t) 

M 2 = If, -I) = - iy) t> + y/l\z I) 

M 3 = \h\) = ^\{x+w)\)-\[i\zV 

u 4 = ||,-|) = 



(3) 



12' 2/ v^l 

were taken as basis functions to deal with the present KL model Hamiltonian j8]. The 
functions \x), \y), \z) have the lattice periodicity, and transform the atomic orb itals p x ,p y ,p z , 
respectively. Here | f), | \), stand for the spin eigenf unctions. The order of (J3j) , follows that 
proposed in Ref. which is: hh + 3,lh_i,lh + i,hh_3. 
We have introduced 

W u = A lK 2 T + V(z) 
W22 = A 2 k 2 t + V(z) 

M = - — (71 + 72) 
M = 7^ — (71 -72), 

with 7j the Liittinger parameters, and mo the bare electron mass. The k t components k XtV , 
are set in-plane respect to the heterostructure interfaces. In Q, J4/O4 stands for the (4 x 4) 
identity/null matrix, while in (|2]), &(z) is an envelope function. Though moderately rough, 
assertion (|2| represents a reliable-accuracy approximation to the V e s, we are interesting in. 
Let us consider a periodic three-layer [A-cladding left (L) layer /B middle (M) layer/ A- 
cladding right (R) layer] heterostructure, in the absence of external fields or strains. In the 
bulk cladding layers, hh and Ih modes mix due to the k p interaction, while the middle slab 
represents a inhospitable medium for holes. At zero valence-band mixing, one has 
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GaP-AlN-GaP, k x = k y = ^f, Ih (a) 
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FIG. 2. Panel (a)/(b) presents the metamorphosis of the effective potential profile V e s for Ih/hh 
(red/blue lines), as a function of k t and layer dimension for a GaP-AlN-GaP heterostructure. 

being 9 a steplike function, and V A / B the potential of the cladding/middle layer. 

Assuming a plane- wave-like dependence on k z for *&(z), and based on the quadratic 
eigenvalue problem (QEP) method, [8] it can be cast 

det [Q(k z )\ = q k 8 z + gi k 6 z + q 2 k\ + q 3 k 2 z + q 4 , (5) 

where Q(k z ) is a second-degree matrix polynomial on the ^-component wavevector k z and 
qi are function of Liittinger semi-empirical valance-band parameters and the components of 
/%.. By dealing with (J5]) , it is straightforward to follow whereas k z is oscillatory or not, and 
thereby the kind of Kff the holes interplay with. 

III. NUMERICAL RESULTS AND DISCUSSION 

On general grounds, for k t the effective band offset is given by the difference for 3D 
band edge levels, [3 E] By letting grow k t , the band mixing effects arise and the effective 
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FIG. 3. Panel (a)/(c) displays a cut of the effective potential profile V^f for Ih/hh (red/blue 
lines), at the interface plane between left and middle layers, as a function of k t . Panel (b)/(d) 
shows the progression of the band offset at the same interface for Ih/hh (red/blue lines), i.e. the 
difference between the upper-edge and lower-edge of the V e Q profile. We have considered a GaP- 
AlN-GaP layered heterostructure. Notice the coincidence of the curves along [10] (solid line) and 
[01] (dashed-dotted line), in-plane directions. 



band offset changes. A comprehensive analysis of this subject, describes larger reduction 
for the piecewise constant effective barrier height V e s as a function of K r , for lh respect to 
that for hh [8]. The mechanism responsible for this behavior is the increment of the term 
n\/m* hh i h (z), yielding even to invert the roles of wells and barriers [6J [Tj as can be derived 
from equation Some authors had declared a shift upward in energy, of the bound states 
in the effective potential well as the transverse wave vector increases, [TT]. Recently a valence- 
band mixing first-principle theory within the EFA was proposed, which approximated the 
superlattice potential energy by considering only the linear and quadratic responses to the 
heterostructure perturbation [12]. 

Today the III-V binary- and ternary-compound semiconductor alloys, continue attract- 
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FIG. 4. Panel (a)/(b) displays the metamorphosis of the effective potential profile V e Q for Ih/hh 
(red/blue lines), as a function of n T and layer dimension for a GaN-AlSb-GaN heterostructure. 

ing interest, for real-life applications and facilities. Nano-dimensional systems based on 
AlAs, AlSb, and AlGaN, together with quantum wells composed on GaAs, InAs, GaP, 
GaSb and GaN components, are a robust platform for material science, condensed-matter 
physics and for a rapidly raising field of Spintronics [13]. In this concern, we briefly describe 
here some phenomenological properties of several representative III-V binary-compound 
semiconductor alloys. To gain some insight into the rather complicated influence of the 
band mixing parameter k t , on the effective band offset, we display several graphics in the 
present section. The central point here, is a reliable numerical simulation for the spatial 
distribution of V e Q while the valence-band mixing increases from k t m (uncoupled holes) 
to k t = O.lA -1 (strong hole band mixing). This purpose requires a solution of ([2J looking 
for a systematic start-point theoretical treatment of highly specialized III-V semiconductor 
binary-compound cases of interest. We have set a width of 25 A for the external cladding- 
layer L and R, while for the middle one we have taken a thickness of 50 A. 

Figure [I] demonstrates that the standard rectangular distribution for V z Q, is a consistent 
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potential-energy trial for hh (blue lines), applicable in the wide range of k t [see panel (b)]. 
On the contrary, panel (a) remarks that the rectangular shape for V z is no longer valid for 
Ih (red lines), as k t increases. It is in this very sense, when the valence-band mixing effects 
get rise, that become unavoidable to refer an effective band offset for a realistic description 
of the interplay of the envisioned physical structure with holes. We display in panel (a), the 
metamorphosis of V e g for Ih, as a function of k t and layer dimension for a GaAs-AlAs-GaAs 
heterostructure. Two changes are neatly observable, namely: the energy edge of both left 
and right cladding-layers steps up in almost 0.5 eV, while for the middle one it remains 
almost constant. The last departs from the V e R evolution for hh, where all borders move 
up almost rigidly [see panel (b)]. Let us, briefly comment these results in an attempt of a 
qualitative comparison with ones given in j3], when a finding of hole states for a quantum 
well embedded in asymptotically flat potential was addressed. We conjecture the obtained 
by them dispersion of lower- hh bound levels [3], to be considered reasonably correct even for 
larger k t , due to the quasi-flat profile of V e g, observed for hh in our Figjljb), at any value 
of band mixing. While for Ih states, some changes could be expected in the lower bound 
states [3], as the effective barrier height — for Ih — , is no longer fixed, but rather becomes a 
dynamically modifiable quantity with k t , as can be seen from our FigjT](a). Although not 
shown here owing to brevity, a similar behavior was found for other middle-layer alloys (AlSb, 
A1P, A1N). Described above features, remain under modification of the in-plane direction. 

Figures [2] and [3] illustrate the phenomenology in another group of III-V semiconductor 
binary-compounds. They show the V e s metamorphosis for a GaP-AlN-GaP layered het- 
erostructure. In panel (a) of Fig|2j we observed a reduction of the effective band offset for Ih 
(red lines), similar to that of a GaAs-AlAs-GaAs heterostructure [see Figjl](a)]. However, 
the panel (b) of Fig J5] illustrates a particularly different behavior for hh (blue lines), due to a 
decrease of the energy edge for both left and right cladding-layers. A more explicit represen- 
tation of that, is shown in Figj3^c). Concerning results of all panels in Fig|3j we found useful 
to underline the coincidence for V e g profile as well as for band offset along [10] (solid line) 
and [01] (dashed-dotted line), in-plane directions. The last is a clear evidence of a partial 
isotropic behavior for the overlapped lines. From Figj3](c), it can be noticed, how V e g dimin- 
ishes its lower-edge profile down to —0.05 eV roughly, along the [11] in-plane direction (blue 
dashed line). Meanwhile, the upper-edge profile evolution plotted in Fig|3^c), goes on nearly 
constant along the same direction. On the other hand, in this kind of heterostructure, the 
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FIG. 5. Root-locus-like map for the eigenvalues k z from QEP's solutions (j5j), as a function 
of kx = k y with k t € [10 -6 ,10 -1 ] A" 1 for: (a) AlAs, (b) AlSb, (c) A1P and (d) AIN binary 
compounds. We had assumed V{z) = 0.498 eV, and the longitudinal energy E = 0.475 eV. The 
map evolves outward coordinate center of imaginary axis, with different values of k z for both Ih 
and hh. 



in-plane full-isotropic performance of V e g evolution, is absence. This fact is observed from 
panels (c)-(d) of Figj3] However, it can be assumed again a rectangular-like distribution of 
V z Q for hh (blue lines), as an acceptable reference frame along the [10] (solid line) and 
[01] (dashed-dotted line), in-plane directions. Moreover, notice that the difference between 
profile edges [see Fig|3](d)] is less than 0.1 eV at the higher n T value, respect to the zero 
hole-mixing point. Concerning the same in-plane directions, worthwhile to underline that 
for Ih, a slightly reduced V c s profile progression, stays as trend [see Figj3|a), red solid line]. 
Notice at panel Figj3](b), that the V e g band offset evolves up to 0.4 eV along [10] and [01] 
in-plane directions at maximum k t , while it reduces almost to 0.2 eV for the [11] interface 
orientation. 

Figure [4] plots the evolution of V e s for a GaN-AlSb-GaN layered heterostructure. In the 
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presence of nitrides in the external layers, the noticeable increment of the effective band 
offset for Ih (red lines) observed in panel (a) of Fig|4j contrasts with the reduction we had 
shown in panel (a) of Fig|3]for non-nitride alloys. Meanwhile in this case, the hh (blue lines) 
feel a similar distribution as in GaAs-AlAs-GaAs heterostructure (see Figure [I]). 

For the sake of illustration, we retrieve a procedure that relies on a robust classical con- 
trol theory technique, a so-called root locus pH]. Once we have quoted the eigenvalues k z 
of (pH), it is then straightforward to generate a plot in the complex plane, symbolizing the 
locations of k z values that rise as a band mixing parameter k t changes. We take advantage of 
the root-locus-like know-how, to promptly identify evanescent modes, keeping in mind that 
complex (or pure imaginary) solutions are displayed in the up(down)-half plane (see Figure 
[5]). These values are forbidden for some layers and represent unstable solutions, underlining 
the lack of hospitality of these slabs for oscillating modes. Although not shown here, the 
opposite examination is also suitable for propagating modes, which become patterned in 
left(right)-half plane, and are equated with stable solutions for given layers. Thus we are 
able "to stamp" on a 2D-map language, the stability domain analysis of the envisioned het- 
erostructure under a quantum-transport problem. This way, we are presenting an unfamiliar 
methodology to deal with solid-state low- dimensional physical phenomenology. 



IV. CONCLUSIONS 

We present an alternative procedure to simulate graphically, the phenomenon of the trans- 
verse degree of freedom dependence of the effective scattering potential. For low-intensity 
valence-band mixing regime, a fixed-height rectangular distribution of the potential-energy, 
is a good trial as a standard reference frame, for a theoretical treatment involving both fla- 
vors of holes in the materials under study. However, this assertion is no longer valid for light 
holes, whenever the mixing effects reveal, and a mutable effective band offset profile has to 
be taken into account, whenever the in-plane kinetic energy changes. Meanwhile, for heavy 
holes the effective band offset profile assumed at k t « remains robust for finite values of 
valence-band mixing. Evidences of this sort, are relevant in experimental applications, and 
in theoretical analysis of coupled-hole tunneling. 
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